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ABSTRACT 

Stick-slip - the sequence of mechanical instabilities through which a slider advances on a solid 
substrate - is pervasive throughout sliding friction, from nano to geological scales. Here we suggest 
that trapped cold ions in an optical lattice can also be of help in understanding stick-slip friction, 
and also the way friction changes when one of the sliders undergoes structural transitions. For that 
scope, we simulated the dynamical properties of a 101-ions chain, driven to slide back and forth by a 
slowly oscillating electric field in an incommensurate periodic "corrugation" potential of increasing 
magnitude Uq. We found the chain sliding to switch, as Uo increases and before the Aubry transition, 
from a smooth-sliding regime with low dissipation to a stick-slip regime with high dissipation. In 
the stick-slip regime the onset of overall sliding is preceded by precursor events consisting of partial 
slips of few ions only, leading to partial depinning of the chain, a nutshell remnant of precursor 
events at the onset of motion also observed in macroscopic sliders. Seeking to identify the possible 
effects on friction of a structural transition, we reduced the trapping potential aspect ratio until the 
ion chain shape turned from linear to zigzag. Dynamic friction was found to rise at the transition, 
reflecting the opening of newer dissipation channels. 



I. INTRODUCTION 

Similarly to colloidal monolayers driven across laser- 
generated surfaces^ ^, linear chains of cold ions trapped 
inside optical lattices have been recently proposed as 
novel candidates for studies in the field of friction^. One 
of the motivations has been the possibility to observe, 
thanks to their exceptional parameter tunability, the long 
theorized Aubry transition, namely the switch between a 
regular frictional state and the "superlubric" state of van- 
ishing static friction between idealized incommensurate 
one dimensional (ID) "crystals". The key feature of fric- 
tion between solid bodies is hysteresis, that is the differ- 
ence between to and fro motion. In time-periodic sliding 
motion for example, hysteresis is responsible for the finite 
area enclosed by the force-displacement cycle, which ex- 
actly equals the frictional heat per cycle. Smallest when 
the sliding regime is smooth, friction turns large when 
sliding occurs by stick-slip - a discontinuous stop and go 
which constitutes the largest and commonest source of 
frictional hysteresis. Generally triggered by mechanical 
instabilities, stick-slip takes place at geological, ordinary, 
and at nanometer length scales alike^ ^. Restricting here 
to the nano and microscale, which is the focus of much 
current work, we are naturally interested in microscop- 
ical systems exhibiting a controlled transition between 
smooth and stick-slip sliding regimes. 

One dimensional periodic ("crystalline") sliding mod- 
els, although highly simplified, have long been used to il- 
lustrate frictional phenomena between periodic latticed. 
In the so-called Frenkel Kontorova (FK) model^, a har- 
monic chain of classical masses with average spacing ao in 
a sinusoidal periodic potential of amplitude Uq and wave- 
length A (leading to a commensuration ratio r]=ao/X be- 



tween the two) idealizes the sliding of two crystalline sur- 
faces. Irrational values of 77 characterize the most inter- 
esting incommensurate case between slider and substrate. 
Aubr}i^^ proved long ago that a transition (where the 
ground state "hull function" exhibits analyticity break- 
ing) occurs for increasing /7o, from what is now known as 
a superlubric state where the static friction Fs - the min- 
imal force capable of initiating sliding - is exactly zero, to 
a pinned state where Fs is finite. While exceptionally low 
friction between incommensurate 3D surfaces has indeed 
been observed experimentally^^, there had not been so 
far experimental demonstrations of the Aubry transition 
in genuinely ID systems. Cold ion traps were recently 
invoked as possible candidates to display the Aubry tran- 
sition, thereby surprisingly entering the field of nanotri- 
bology^. Although not identical to the FK model, the 
physics of repulsive ID particles is expected to be es- 
sentially the same as each particle can still be seen as 
occupying the center of some overall harmonic potential. 
Experimentall}'^^, chains of up to several tens of positive 
ions such as Ca+ can be stabilized using rf quadrupo- 
lar fields and cooled down to temperatures below 1 jaK. 
By tuning the confining cigar-shaped potential to a suffi- 
ciently elongated form, the ions can be forced to form lin- 
ear chains. The periodic optical lattice potential for the 
ions is provided by a laser standing wave (see Fig. [l^). 
The confined ion chains do constitute ID crystal seg- 
ments, but are not really homogeneous. The nearest 
neighbor ion-ion distance, fairly constant at the center, 
increases at the periphery and diverges near the extrem- 
ities, as shown in Fig. [TJd. Still, the chain center is a 
reasonable realization of an FK-like model, and some 
of the properties of an ideal infinite system can be in 
principle realized and observed there. When A is incom- 
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FIG. 2. (Color online) Effective static friction of a 101- 
ions chain in the linear (red) and zigzag (blue) configuration 
plotted against the corrugation potential amplitude Uq. The 
black dashed line is the static friction of the ideal FK model. 
(From Ref. g]) 



FIG. 1. (a) Schematic of a linear ion chain trapped by an 
anisotropic confining potential. Typical ions used in experi- 
ments are Ca^ or Mg^, the ion- ion distance at the center of 
the chain is of the order of few /xm. (b) Nearest-neighbor dis- 
tance (in dimensionless units, see Section [ll]) between the ions 
of a 101-ions chain at rest and in absence of the corrugation 
potential, (c) A possible experimental set up for the study of 
the dynamics of the chain using an external oscillating electric 
field. 



mensurate with respect to the central ion-ion spacing 
one can achieve, according to our recent predictions^, 
a strong and observable remnant of the Aubry transi- 
tion also in such trapped ion chains. In the confined ion 
chain, the standard Aubry transition, which in the infi- 
nite chain occurs when the periodic potential ( "corruga- 
tion") amplitude Uq exceeds some critical threshold Uc^ 
is replaced by a static, symmetry breaking transition of 
the ground state chain configuration and geometry. Be- 
nassi and coworkers^ proposed to observe this transition 
by measuring the external uniform force Fn needed to re- 
store the symmetry. Simulations indeed showed that the 
effective static friction force Fn behaves and grows as a 
function of Uq > Uc very closely like the static friction 
force Fs of the ideal infinite chain thus demonstrating 
the connection between the two (see Fig. [2| . 

In this work we move from static friction to the dy- 
namical sliding properties of the ion chain, once depinned 
by an additional external electric field, as a function of 
the periodic corrugation amplitude Uq (see Fig. [l]^). We 
show that the trapped cold ions can slide either smoothly 
or by stick-slip, with a parameter-controlled transition 
and a correspondingly strong frictional rise between the 
former and the latter. Drawing an analogy with macro- 
scopic frictional experiments, the corrugation Uq plays 



here the role of the load in ordinary sliding friction. As 
expected, we find and characterize the transition from a 
poorly dissipative smooth sliding regime to a highly dis- 
sipative stick-slip regime as Uq is increased. Moreover, 
since the cold ion ground state geometric configuration 
can be pushed across parameter-driven structural transi- 
tions by changing the trapping potential conditions, we 
investigate the effect of a phase transition on sliding fric- 
tion, which is of interest as welP^. As is known both 
theoreticall}ff^lI^ and experimentally^, a change of as- 
pect ratio in the confining trap effective potential causes 
the ion chain to cross a series of structural transitions. 
If and when for a sufficiently long chain these transitions 
can be considered continuous, the friction behavior near 
the transition point could show remnants of the chain's 
critical behavior, as recently suggested theoretically^. 



In anticipation of future experiments, we carried out 
classical molecular dynamics simulations of a 101-ions 
chain sliding in a golden ratio incommensurate corru- 
gated potential, with a view to predict and discuss the 
basic dynamic frictional phenomena of an electric field so- 
licited trapped ion chain. In Section [TT| we will describe 
the model and the protocol used for the MD simulations. 
Section |III| will be devoted to the resulting smooth to 
stick-slip frictional switch and the observation of precur- 
sor events at the onset of the chain sliding. In Section IV 
we describe the change of the frictional behavior across 
the linear-zigzag structural transition. Finally Section [V| 
contains our discussion and conclusions. 
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FIG. 3. 
chain, r, 

z axis, (b) Zigzag-helix transition. The order parameter is 
the number of windings of the ions around the z axis. 



II. TRAPPED ION CHAIN MODEL AND 
SLIDING SIMULATION PROTOCOL 

The effective potential of an ion of charge in a linear 
anisotropic (Paul) trap can be written a^^ 
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where m is the mass of the ion and uj\ and cjy are the 
strengths of the confining effective potential, supposed 
to be harmonic, in the transverse and longitudinal direc- 
tions. In order to work in dimensionless units we define 
the length unit d\ 
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We then measure masses in units of m, charges in units 
of energy in units of g^/(47reo<i), forces in units of 
g'^/(47reo(i^) and time in units of 1 / uj±. The effective 
Hamiltonian of N trapped ions is therP^^ 

Heff = Eif + ^[c.i(x,2 + yf) + o.jzf] (3) 

i=l 

where the sinusoidal term represents a laser induced peri- 
odic potential, mimicking the corrugation of a hypothet- 
ical crystalline substrate lattice. 

The ground state geometry of the ions at T=0 depends 
on the aspect ratio R={uj\\/uj±)'^ of the anisotropic har- 
monic confining potential. For small enough R the po- 
tential is cigar-like, and ions form a linear chain along 
the trap symmetry axis, z. As R is increased there is a 
sequence of shape transitions: first from a straight chain 
to a planar zigzag chain; next, a second transition where 
planarity is lost, and the planar zigzag turns into a he- 
lix. Still at T=0, and for an infinite chain, both classical 
transitions are continuou^^ as shown in Fig. [sj 

We simulated, following Benassi et al.^, a chain of 
101 positive ions choosing cJm =0.0005, i?=0.0005 and 



r]=X/ao=2/{l -h \/5), where is the center ion- ion spac- 
ing. Chain sliding is caused by an external slowly os- 
cillating electric field 'E{t)=zEosin{Qt) acting on each 
ion in the longitudinal direction z. We carried out clas- 
sical damped molecular dynamics (MD) integrating the 
equations of motion using a standard velocity- Verlet al- 
gorithm with time step At=0.005. At each time the total 
force acting on the i^^ ion is given by: 
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where we have respectively the force due to the ion-ion 
Coulomb repulsion, the confining potential force, the cor- 
rugation force and a velocity dependent dissipative force 
controlled by a damping parameter 7. There is no ran- 
dom force, corresponding to our "T=0" background as- 
sumption. The trap confinement plus oscillating poten- 
tial is given by 



zEosin{Qt)^ 



(5) 



which is a confining parabola of ver- 
tex Ztrap=Eosin{Qt)/uj'^^ moving at velocity 

Vtrap=Eo^cos{rtt) /u'^^. In ordcr to follow stick-slip, 
when present, we monitored the distance of the center 
of mass of the chain from the minimum of the moving 
parabola: 



^^cm(^) ^trap(^) ^cm(^) 
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Fi gure [4] shows an example of the time evolution of Szcfn 
corresponding to a sequence of external electric field oscil- 
lations. The dynamic friction of the system is computed 
as the work W done by the oscillating electric field on all 
the particles: 



1 ^ r 

Wk = j^^ dt[^^ • zEosin{Qt)], 

i=i 



(7) 



Where the integral is calculated over the k^^ period cor- 
responding to the electric field going from its minimum 
value —Eq to its maximum value +^0 and back. The 
final estimate of the dynamic friction is obtained from 
the average of the M samples Wk measured during the 
whole trajectory: W=j^ 



III. RESULTS: DYNAMIC FRICTION, 
STICK-SLIP AND PRECURSOR EVENTS 

Figure [5^ shows the frictional work W done by the 
external electric field on the trapped ion chain as a func- 
tion of the corrugation amplitude Uq. The inset depicts 
schematically the expected behavior in the infinite in- 
commensurate FK model, where for Uq < Uc^ below the 
Aubry transition, motion takes place without static fric- 
tion, and kinetic friction vanishes in the limit of infinitely 
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FIG. 4. Dynamics of 6zcm (see text) under periodic sliding. 
At t=0 the chain is almost symmetric with respect to the 
vertex of the confining parabola and 6zcm ~ 0; if the chain is 
locked to the substrate Szcm increases as the confining poten- 
tial moves until a slip event occurs, which corresponds to a 
sudden drop to zero of 6zcm- The plot shows a result from a 
simulation where many slip events are observed during each 
oscillation back and forth of the chain over the substrate. The 
grey dashed line corresponds to the external electric field. 



slow sliding. In the chain of trapped ions, finite and in- 
homogeneous, the static friction force needed for overall 
chain motion is nonzero for all corrugations, since the 
two extremities are always locked to the corrugation po- 
tential. Correspondingly, there is upon sliding a finite 
frictional dissipation W for all values of Uq. For the 
chosen external field frequency 1^=0.0002 and amplitude 
£^0=0.1625 friction grows steadily with corrugation Uq 
and nothing significant happens to the dynamic friction 
W across the nominal^ Aubry value Uc=0.0W2S. In fact 
the dynamics of the chain does not change appreciably 
above Uq ~ 0.01. 

We can resolve, based on the detailed nature of the slid- 
ing trajectories, three different dynamical regimes. Fig- 
ure (5)3 shows the dynamics of Szcm for three representa- 
tive values of Uq. For small corrugations Uq ^ /7oi =0.005 
the chain follows smoothly the external force and the fric- 
tion is modest (with a value determined by, and grow- 
ing with, the sliding velocity, in turn proportional to 
Q). As Uq is increased further, the smooth sliding dy- 
namics is replaced by a regular, time-periodic stick-slip 
regime with accompanying increase of dissipation. So 
long as Uqi ^ Uq ^ /7o2=0.01 the slip magnitude is fairly 
constant during each oscillation. For larger corrugations 
finally, Uq ^ U02 the chain enters a chaotic regime of 
irregular stick-slip. 

Figure [6] shows details of the trajectories of the chain 
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FIG. 5. (a) Dynamic friction W of the oscillating ion chain 
as a function of the corrugation amplitude Uq. The values 
of ^0=0.1625, ^^=0.0002, 7=0.01 were chosen so as to yield 
4-5 slip events during each oscillation in correspondence of 
the highest values of Uq investigated, (b) Dynamics of Szcm 
in the three different regimes of smooth-sliding, regular and 
irregular stick-slip. The time interval plotted corresponds ap- 
proximately to half a period of the oscillating electric field 
which increases from —Eq to -\-Eo. 



center of mass and of all individual ions in the three 
regimes. For Uq ^ Uqi the chain is weakly pinned at the 
inversions of motion occurring for E(t)=±£^o- After the 
depinning the chain follows smoothly the external force; 
small oscillations of Szcm are due to internal motion of the 
chain (see Fig. |6^,b). A number of interesting features 
appear at the onset of stick-slip Uq ^ Uqi and they are 
shown in Fig. [6]:-f. The head and the tail of the ion crys- 
tal are locked to the corrugation potential (see top and 
bottom parts of Fig. |6]i, initial times) while the truly in- 
commensurate central part is free to slide, thus increasing 
the ion density of the head, reducing that of the tail and 
inducing partial depinnings of the chain. The precursor 
events appear in the head part (see top part of Fig. [6]i, 
times between t=1.16 and t=1.22). Although different 
and connected with inhomogeneity of stress rather than 
of contact, partial precursors were also shown to pr ecede 
the onset of macroscopic sliding by Fineberg's groupP^I^. 
Partial slips of the chain always start within the central 
superlubric region and proceed moving in the direction of 
the external force. They are present only in the stick-slip 
regime and disappear by increasing the average pulling 
velocity \vtrap\='^EoQ/7ruj?; of the trapping potential (see 
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Fig. [7|. Following the precursors, as stress accumulates 
more and more in time, the system undergoes a mechan- 
ical instability typical of stick-slipP. The ensuing main 
slip of the whole system is triggered by the creation of a 
kink-antikink pair (see Fig. [8|. In the regular stick-slip 
regime Uqi ^ Uq ^ U02 the tail drives the sliding onset, 
and the triggering pair forms in the tail region, while the 
chain center and front are still free to slide (see Fig. |6]i, 
times between t=1.22 and t=1.25). In the chaotic stick- 
slip regime Uq ^ U02 instead, the onset of sliding is dif- 
ferent. The central superlubric flow and partial slips of 
the front ions first bring the chain into a metastable state 
where each ion is locked to the corrugation (see Fig. [6]p, 
times between t=3.5 and t=3.6). The generation of a 
kink-antikink pair forming now in the chain head, as op- 
posed to the chain tail of the previous regime, eventually 
leads to global sliding, as shown by the main slip event 
at the initial times of Fig. [Bf. This tail-to-head switch 
of the triggering event is a characteristic signature al- 
ways accompanying the passage from regular to chaotic 
stick-slip. On the other hand, neither the sliding onset 
dynamics nor the dynamical friction magnitude finally 
display any particular feature or singularity when the cor- 
rugation grows across the Aubry transition. This result 
underlines a substantial difference between the frictional 
behavior of this short and inhomogeneous chain, and that 
expected of an ideally infinite and uniform FK-like chain. 
In the latter and ideal system there is no other transition 
than Aubry, and in the limit of zero sliding speed stick- 
slip sliding can only take place when static friction turns 
nonzero, which is above the Aubry transition. Figure |9] 
shows details of the sliding dynamics in the strongly cor- 
rugated, chaotic stick-slip regime. The average nearest 
neighbor distance between the central 31-ions portion is 
displayed as a function of time between two main slips of 
the crystal. The initial passage of the "superlubric" front 
brings the central part to a commensurate configuration 
with ions spaced exactly by 2 A (instead of the original 
golden ratio spacing 1.618 A) from one another. This 
dynamically induced commensuration brings the whole 
chain to a temporarily locked state thus increasing the 
static friction force needed for the onset of overall mo- 
tioiP^. Subsequent depinning of the chain off this locked 
state only occurs as the external force grows further, and 
is sudden. This two-stage nature of sliding, and the rel- 
ative abruptness of the depinning is at the origin of the 
chaotic behaviour of stick-slip in this regime. 



IV. FRICTION SINGULARITY AT THE 
LINEAR-ZIGZAG STRUCTURAL TRANSITION 

Benassi et al. recently proposed that interesting fric- 
tional changes or anomalies could be observed in pres- 
ence of structural phase transition^. Although a 101- 
ions chain is a long way away from an infinite system, 
it is still interesting to find out what singularities would 
friction develop upon an overall, collective shape change. 



To study that, we carried out simulations at values of 
the trapping potential aspect ratio R straddling the crit- 
ical value Rc=O.OOOS for the linear-zigzag transition. Al- 
though this kind of effect should be quite general, the 
expected delicacy of this frictional feature should best 
become apparent under sliding conditions with limited 
noise, such as those expected at weak corrugations. Set- 
ting /7o=0-0008 and also 7, Q and Eq values which lead 
to a smooth and gentle sliding dynamics (disturbing the 
chain to a minimal extent) we obtained the dynamic fric- 
tion of Figure [10^. Deep enough in the linear chain 
regime {R < 0.00064) the ions remains in a strictly ID 
configuration during the whole dynamics. Here, only 
longitudinal internal vibration degrees of freedom of the 
chain are excited and W is small and essentially indepen- 
dent of R. As the critical anisotropy Rc is approached 
the transverse vibration modes of the chain soften and be- 
come rather suddenly excitable, a new dissipative channel 
opens and W rises anticipating the linear-zigzag transi- 
tion. In a hypothetical infinite chain, where the transi- 
tion occurs continuously and critically, the frictional be- 
havior will presumably also exhibit a critical singularity. 
Given the finite chain size, the frictional rise is smooth, 
although it can still be sharpened by reducing the cor- 
rugation amplitude. Figures [To]3,c show results of sim- 
ulations with Uo reduced down to 0.0004, 0.0001 and a 
lower value of Eq. The transverse mode excitation onset 
occurs nearer to Rc. 

That result is made clearer in Fig. [lO]i where the stan- 
dard deviation a^^^ of the maximum displacement of the 
ions away from the z axis is plotted against the aspect 
ratio R. Denoting with < . . . > the time average over 
the whole trajectory, ar^y and rxy are defined as: 

^r^y = \J< i^xy- < r^y >)2 > (8) 

r^y = Max|i^i,Ar_}(^xf +^2) (9) 

(Jr^y is zero until the transverse modes begin to be excited 
and it displays a maximum near the critical point Rc. 

A further increase of R drives a subsequent transition 
of the trapped ion geometry from zigzag to helix. Fig- 
ure [11] shows the behavior which we obtained for the 
dynamic friction across this transition. Not much hap- 
pens here. In this case in fact, on both side of the critical 
point the chain is already in a 3D configuration and its 
dissipative properties are not significantly affected by the 
weaker helical distortion. 



V. CONCLUSIONS 

By means of classical damped MD we simulated a 101- 
ions linear chain executing a forced time-periodic sliding 
over a (laser induced) space-periodic "corrugation" po- 
tential of strength Uq whose wavelength A was golden 
mean incommensurate with respect to the center ion-ion 
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FIG. 6. (a) Dynamics of bzcm in the smooth sUding regime; panel (b) shows the corresponding trajectories of each ion. (c) 
Dynamics of bzcm in the regular stick-slip regime; panel (d) shows the corresponding trajectories of each ion. (e) Dynamics of 
bzcm in the irregular stick-slip regime; panel (f) shows the corresponding trajectories of each ion. 
The features highlighted in the figure are discussed in Section [III] 



spacing. As [/q was increased the system turned from 
smooth sliding to a stick-slip sliding regime, first regular 
and then chaotic. We observed as expected an increase 
in the dissipation rate at the onset of stick-slip, as is also 
observed in macroscopic dry friction when the loading 
force is increased. 

Three separate frictional regimes were identified, as 
a function of corrugation amplitude. A smooth sliding 
one for weak corrugation was followed by time-periodic 
stick-slip sliding at larger corrugation, eventually leading 
to chaotic stick-slip for even larger corrugation. 

Due to the inhomogeneity of the ion crystal the fric- 
tional dynamics of the ion chain showed several novel fea- 
tures reflected in the ion trajectories. The two chain ex- 
tremities were always pinned while the incommensurate 
central part was free to slide following the external force. 
The onset of motion in the stick-slip regime was charac- 
terized by the presence of precursor events, i.e. partial 
slips of side chain portions induced by the superlubric 
flow of the truly incommensurate central part. 

The chaotic stick-slip at large corrugation was con- 
nected by an interesting two-stage process. First, super- 
lubric sliding of the central portion brought the chain to 
a temporary commensurate state, locked to the periodic 
corrugation. Subsequently, as force grew, a kink-antikink 
pair was generated and propagated toward the extremi- 
ties eventually inducing the slip of the whole system. 

We also studied the possible anomalies of friction dis- 
sipation across the structural phase transformations of 
the trapped ions obtained by varying the aspect ratio 



R={uj\\juj\})^ of the harmonic trapping potential. At 
fixed c^ii , as R was increased, the ion chain transformed 
first from a linear configuration to a planar zigzag and 
then to a helix. The energy dissipation increased charac- 
teristically at the linear to zigzag transition. Conversely, 
the zigzag to helix transition did not yield significant fric- 
tional changes. 

Considerations about experimentally accessible system 
parameters, which are given in the appendix, suggest that 
some of these features, if not all, should become observ- 
able in future experiments with Ca+ ion traps. 
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Appendix: Practical Parameter Choices 

In order to create an optical lattice the laser wave- 
length must fit one of the electronic transition of the 
chosen ions. For Ca+ the S1/2-P1/2 transition at 397 
nm is naturally exploited, leading to a lattice spacing 



7 



Uo=18x10"^ |Vtrap|=0.01 Y=0.015 



Uo=6x10"^ |Vtrap|=0.04 Y=0.01 



N-400 




36 37 38 39 40 41 42 43 44 

tx10^ 

Uo=18x10-2 |Vtrapl=0-1 Y=0.015 



N-200 




1.65 1.80 1.95 2.10 2.25 2.40 
x3 



tx10' 



FIG. 7. Precursor events observed between two main slip 
events for (a) \vtrap\='^0~^ and (b) \vtrap\='^0~'^ . As the 
pulling velocity is increased the precursor events disappear 
signaling the transition from stick-slip motion to smooth slid- 
ing. 



A ~ 200 nm. An achievable amplitude for the corruga- 
tion potential is of the order of 10~^^ J. If we consider a 
transverse trapping frequency uj±=27tx4: MHz we get (in 
dimensionless units) A=0.115, Uo=2.31x 10~^. A practi- 
cally achieved temperature is T=10~^ corresponding to 
1 fiK. Setting an aspect ratio i?=0.0005 we have ~16 
A, therefore ions are separated by several lattice spac- 
ings, corresponding to a much "weaker" kind of incom- 
mensurability that the golden ratio used in the study so 
far. Moreover for such a small value of Uq the chain is 
almost free to slide, in this case also the two extremi- 
ties being weakly anchored to the substrate. Stick-slip 
is therefore expected to occur only at very small pulling 
velocities when, after each slip, the chain is allowed to 
relax in a new metastable pinned configuration. We per- 
formed simulations at the experimental parameters for 
the 101-ions and a 35-ions chain using the same proto- 
col described above. We chose 7=0.0005 and we used 
a Langevin thermostat for the simulations at finite T. 
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FIG. 8. Example of a triggering event inducing the main slip. 
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FIG. 9. Average ion- ion distance (in units of A) versus time; 
we considered only the central 31 ions of the chain. After the 
passage of the superlubric front the center of the chain is left 
in a commensurate configuration with ions equally spaced by 
2 A. The same kind of dynamics occurs at the onset of motion 
at each values of Uo in the stick-slip region. 



We also set Q and Eq in order to test different average 
velocities of the moving confining parabola. 

Let us consider first the 101-ions chain. At very low 
pulling velocity we observed stick-slip motion, the slip 
amplitude being of the order of a few lattice parameters 
(see Fig. p!2^). The dynamics of the slip events is sim- 
ple and no longer shows the variety of features described 
previously for the golden ratio incommensurability. The 
pinned chain remained stable upon pulling until a weak 
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FIG. 10. (a) Dynamic friction across the linear-zigzag tran- 
sition for t/o=0.0008. Parameters used are: 7=0.0005, 
£;o=0.005078125, Q=0.0032. Without corrugation potential 
the chain moves as a rigid body at a very low dissipation 
rate independent on R (not shown). As Uo is switched on, 
the internal degrees of freedom of the chain begin to dissi- 
pate. In this set of simulations the transverse modes begin to 
be excited at ~ 0.00064. (b) Dynamic friction across the 
linear-zigzag transition for [/o=0.0004 and (c) [/o=0.0001, us- 
ing ^0=0.002539063, Q=0.0032, 7=0.0005. By lowering the 
external electric field and the amplitude of the substrate po- 
tential the onset of the excitation of the transverse modes is 
shifted towards the static critical value i?c=0.0008. In this 
set of simulations the transverse modes begin to be excited at 
R ^ 0.00072 and R ^ 0.00078 respectively for [/o=0.0004 and 
[/o =0.0001. (d) Standard deviation ar^y of the maximum dis- 
placement rxy of the ions from the z axis plotted as a function 
of R and measured using the whole trajectories of the simu- 
lations of panels (a),(b),(c). ar^y is zero until the transverse 
modes begin to be excited and it shows a clear maximum near 
Rc, which becomes sharper as Uo is decreased. 



FIG. 11. Dynamic friction across the zigzag-helix transi- 
tion. The dissipative properties of the chain are not affected 
by the helical distortion. Parameters used are: 7=0.0005, 
^0=0.040625, 11=0.0004, [/o=0.0008. 
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compression generated within the chain propagates to- 
ward the head. Simple precursor events are observed 
consisting of partial slips of a central portion of the chain 
at fixed extremities (not shown). Figure 12 3 shows the 
temperature Tcm measured in the center of mass frame, 
displaying the expected inverse sawtooth behavior when 
plotted in a semilogarithmic scale. After each slip Tcm 
increases, and is then exponentially damped by the ther- 
mostat. As the pulling velocity is increased the dynamics 
turns gradually into a smooth-sliding regime, as shown 
in Fig. [I2]3,d. 

Chains of a few tens of ions may be more easily stabilized 
inside a trap. We performed simulations with a 35-ions 
chain, using the same parameters as above. In this case 
the central ion-ion spacing is larger and corresponds to 
ao ~30 A. As shown in Fig. [13^ we observed that stick- 
slip motion is again preserved in this case as well for small 
enough average pulling velocities. The slip amplitude is 
larger than that observed for the 101-ions chain indicat- 



FIG. 12. 101-ions chain. 

(a) 6zcm in the simulation at \vtrap\='2xl0~^ and T=10~'^ 
showing stick-slip and precursor events (indicated by the ar- 
rows) corresponding to partial slips of the central portion of 
the chain, (b) Temperature measured in the center of mass 
frame in the simulation of panel (a). After each slip event 
the internal temperature is raised and then it is exponentially 
damped by the thermostat. Precursor events give rise to the 
smaller peaks indicated by the arrows. This trend is always 
observed in the stick-slip regime, (c) For \vtrap\='^0~^ and 
T=10~^ the dynamics is chaotic. Slips of different magnitude 
occur as the pulling velocity changes during each oscillation, 
(d) Further increasing |'ytrap|=5x 10~^ the dynamics finally 
turns into a smooth-sliding regime, even at T=0. 



ing a stronger pinning to the substrate, as is reasonable 
to expect given a larger prevalence of extremities. From 
the plot of the temperature of Fig. [l3b we see that small 
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adjustments of the inner ions occur, which give rise to 
small peaks in Tcm, however without proper precursor 
events due to the reduced size of the chain. 
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FIG. 13. 35-ions chain. 

{a)6zcm in the simulation at | Vtrap I — 2 X 10"^ and T=10"^ 
An irregular stick-slip regime is observed, (b) Temperature 
measured in the center of mass frame in the simulation of 
panel (a). The highest peaks correspond to the main slip 
events, slip of few ions within the chain give rise to the smaller 
ones, (c) Trajectory of each ions during some slip events of the 
simulation of panels (a) and (b). The chain slip as a whole 
and precursor events are not observed due to the small size 
of the system. {d)5zc7n in the simulation at |i;trap|=2x 10~^ 
and T=10~'^. At the inversions of motion the chain is pinned; 
owing to the small initial pulling velocity a sharp slip event 
occurs after which smooth sliding begins. 
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